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Using the well-known Olami-Feder-Christensen model as our paradigm, we show how to modify 
uniform driven self-organized critical models to generate 1/ f a noise. Our model can reproduce all 
the main features of 1/ f a noise: (1) a is close to one and does not depend on the dimension of the 
system. (2) The l/f a behavior is found for very low frequencies. (3) The spatial correlations do not 
obey a power law. That proves that spatially extended systems based on activation-deactivation 
processes do not have to be point-driven to produce 1//" noise. The essential ingredient is a local 
memory of the activation-deactivation process. 
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A time signal X(t) with zero mean is called 1// Q noise 
or 1/ f a signal if its power spectrum S(f) is proportional 
to 1// Q at low frequencies / with a ss 1. Here, the 
power spectrum is defined as the amplitude squared of 
the Fourier transformed signal, i.e., 
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According to the Wiener-Khinchin theorem, S(f) is the 
Fourier transform of the autocorrelation function C(t) 
which is defined as 
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Consequently, it follows for the autocorrelation function 
of a signal with S(f) oc 1// Q and < a < 1 that 
C(t) oc |t| q . Hence, a 1// Q signal with a close to but 
smaller than 1 is related to (statistical) long-time corre- 
lations which is the reason why 1/ f a noise is considered 
to be a particularly interesting phenomenon a priori. 
Moreover, the omnipresence of 1//" noise in nature is 
one of the oldest puzzles in contemporary physics. It ap- 
pears in a variety of systems from physics, geophysics, 
astrophysics, technology, sociology, and biology. Specific 
examples are the flow of the river Nile jj], sunspot ac- 
tivity pressure variations in the air caused by mu- 
sic and speech || , human coordination Q , and neuronal 
spike trains || . One of the most famous examples is the 
measurement of the voltage drop V on a resistor of resis- 
tance R through which a current I is flowing. The power 
spectrum of the fluctuations around the expected value 
V = RI clearly shows a 1// behavior over many decades 
in the frequency domain [||. 

It is natural to expect that there might be a general 
principle which explains the occurence of 1/ f a signals in 
many of these different systems. However, no generally 
accepted explanation of the ubiquity of 1/f" noise has 
been proposed yet. Indeed, it is possible to find in the 



literature some ad hoc formulas and theories, but most of 
them are based on unverified assumptions, or they catch 
a glimpse of the physics only of some particular system, 
therefore missing to explain the widespread occurence of 
the phenomenon (see Refs. HQ and references therein). 

In 1987 Bak, Tang, and Wiesenfeld (BTW) intro- 
duced the notion of self-organized criticality (SOC) to 
explain the universality of 1// Q noise ||. SOC systems 
are nonequilibrium systems driven by their own dynam- 
ics to a — in a statistical sense — stable state (self- 
organization). Fluctuations around this state, so-called 
avalanches, are characterized by power-law distributions 
in time and space (criticality) implying long-range corre- 
lations (for a recent review on SOC see Refs. |9j]l0|]). This 
automatically leads to a power spectral density exhibit- 
ing a 1//° decay. However, this approach has several 
shortcomings: First and most important, there is no ev- 
idence for power-law s pac e correlations in most systems 
exhibiting l/f a noise |tllfl. This already means that the 
notion of self-organized criticality and l/f a noise is mu- 
tually exclusive in most cases. Second, a is seldom close 
to one in SOC systems. Moreover, the exponent strongly 
depends on the dimension of the SOC system, at least be- 
low the upper critical dimension. Finally, as we will show 
later on, the 1//" behavior in SOC models is observed for 
high frequencies rather than in the low-frequency range. 

Recently, several authors have sucessfully modified 
originally self-organized critical models to overcome these 
problems |l^-|l5(|. Despite the diversity of introduced 
modifications (continuous driving |12| , dissipation [fl3f , 
(quasi-) one-dimensional geometry ]14| , p^[ ), there is one 
common denominator. All these models have a preferred 
propagation direction of the avalanches. This is implic- 
itly defined via specific driving mechanisms. The sys- 
tems in Refs. pl|JTq] are essentially point driven and the 
system in Ref. |12|] is boundary driven. Without these 
special driving mechanisms the models are not able to 
generate 1// Q noise. 

In this paper, we will show that such a preferred prop- 
agation direction is not a necessary condition to obtain 
l/f a noise from SOC models. We propose a simple 
model with uniform driving able to reproduce the above 
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mentioned characteristics of 1/ f a noise. 

One of the main features believed to be relevant for 
the description of 1//" noise is an activation-deactivation 
process jig . This is, for example, realized in stick-slip 
models and, hence, in a model introduced by Olami, 
Feder and Christensen (OFC) in 1992 which was intended 
to mimic the dynamics of earthquakes [ fl7| . In this model, 
a real variable Fi, called stress, is attached to each point 
i of a (/-dimensional cubic lattice of size N = L d . In the 
initial state, the values of F are randomly distributed in 
[0,1] obeying a uniform distribution. The dynamic evolu- 
tion is characterized by slow driving and fast relaxation. 
All sites i = 1, . . . , N are driven at the same rate v as 
long as Fi < F c , i.e. 

F[ = v. (3) 

As soon as one of the Fi 's exceeds the critical threshold 
value F c the stress Fi is redistributed to the 2d nearest 
neighbors of site i, 

F = 0, (4) 

F nn = F nn +0F i . (5) 

Here, (3 describes the level of dissipation. The model is 
conservative for (3 = l/2d and dissipative for < (3 < 
l/2d. The local relaxation continues until all i^'s are 
subcritical again. The sequence of discharges triggered 
in this way is called an avalanche. If more than one site 
is supercritical at any time, the discharges are assumed 
to happen simultaneously. After the avalanche is over 
the slow driving [Eq. (jj|)] sets in again. It is important 
to note that this time scale separation formally implies 
v -► 0. 

In two dimensions, the OFC model is considered to be 
a SOC model provided that open boundary conditions 
are applied. However, it is not clear whether this is only 
true in the conservative case. Recent investigations seem 
to imply that in two dimensions the OFC model could 
only be classified as "almost critical" for values of (3 close 
to but smaller than 0.25 [Q. This is in contradiction to 
claims by other groups that the model is self-organized 
critical even in a certain range of dissipative values of (3 
Jj^ , pj|]2"o| . However, the distribution of avalanches with 
respect to their size obeys a power law in the numerically 
accessable range of system sizes even for a small amount 
of dissipation. 

The model, as it stands, is not a good candidate to 
describe l/f a noise as follows from Fig. [l| There the 
power spectrum of the OFC model is shown for different 
values of (3. The quantity we use as our time signal is the 
avalanche signal 

*(*)=X>*(*-*i), (6) 
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where gj denotes the size (i.e., the number of topplings) 
of the jth avalanche and tj its time of occurence on the 
time scale of the driving. The explicit definition of a 
time scale leads to a time signal with varying temporal 
distances between avalanches. This is in contrast to other 
signals considered so far in the context of SOC and l// a 
noise. The 5 function is motivated by the time scale sepa- 
ration, i.e., we are not able to observe events on the time 
scale of the relaxation. This is also reasonable because 
we are only interested in the low-frequency range where 
1// Q noise is usually found. 

The time signal was recorded after the system had 
reached a stationary state as described in Ref. pj[ . For 
nonzero dissipation, a characteristic frequency occurs in 
the spectrum as already discussed in Refs. Jl^,|2^]. Above 
and especially below the characteristic frequency, there 
is clearly no sign of 1// Q noise. Rather a white noise 
behavior can be identified. 
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FIG. 1. log 10 -log 10 plot of the power spectrum of X(t) in 
the two-dimensional OFC model with open boundaries for 
different fi's and N = 2500, F c = l,v = 0.1. There is clearly 
no sign of l/f a noise. 

This is not in contradiction to the observations in Refs. 
gl]|2|. In Ref. glj, a 1/f behavior was found for fre- 
quencies larger than the characteristic frequency. How- 
ever, a different time signal was used, namely a stress 
signal which is the stress averaged over the lattice sites 
as a function of time. In Ref. p2j, a l//-type behav- 
ior was described. The authors measured the avalanche 
signal in terms of the time scale of the relaxation, i.e., 
the time between different avalanches was essentially set 
to zero. Hence, they observed a high-frequency phe- 
nomenon characterizing the internal temporal develop- 
ment of the avalanches. Therefore, their findings cannot 
be considered as l/f a noise. As we will see later, our 
modifications lead to a 1/f decay in the power spectrum 
below the characteristic frequency which is the range 
where 1// Q noise should be looked for. 

Our model is basically an extension of the OFC model. 
We just add a single new element: The threshold value F c 
becomes a function of space and time mimicking a local 
memory such that each site remembers its (cummulative) 
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history of discharges. The simplest way to model such a 
memory is to implement it by a random process. After 
each toppling the respective F c (i) evolves according to a 
random walk with Gaussian step length 

F c , T {i) = F c , T -i(i) + tnj{i,T), (7) 

where r denotes the number of topplings of site i 
and {r](i,T)} the sequence of uncorrelated normally dis- 
tributed random variables with zero expectation and unit 
variance. The strength of the white noise source is given 
by a. To omit negative threshold values and to con- 
fine the random walk, we impose reflecting boundaries 
at and at F u = 2F C . As initial condition, we use 
F c , (i) = F c . 

Computing the power spectrum for our model, we find 
a clear l/f a decay over several decades for dissipative 
/3's both in one and two dimensions (see Figs. || an d 
|J). The exponent decreases slightly with decreasing f3 
and lies between 0.9 and 1.2. As a rule, the range of the 
1// Q behavior also decreases with decreasing dissipation 
shrinking to zero in the conservative limit. In this limit, 
the crossover leads to a white-noise type of behavior at 
low frequencies. 
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FIG. 2. log 10 -log 10 plot of the power spectrum of X(t) 
in our model with open boundaries for different /3's and 
a = 0.04, d = 2,N = 2500, F c = 1, v = 0.1. The solid line 
with exponent 1.0 is drawn for reference. The dissipative ver- 
sion of the model clearly generates 1/ f a noise. For d — 1 as 
well as for larger system sizes, we obtain similar results. 

Our modification also destroys the power-law distribu- 
tion of the avalanches in the dissipative case. We find an 
exponential distribution instead. Hence, avalanches can- 
not establish long-range correlations anymore through 
the system. This is expected since the responsible mech- 
anism for the critical or "almost critical" behavior is 
marginal synchronization [ p^| , pO[ . This synchronization 
is necessarily destroyed as soon as the threshold varies 
locally (for quenched random thresholds see Ref. 0]). 

Extensive numerical simulations show that our results 
are very stable with respect to variations in the parame- 
ters. Different values of a and F u lead to the same results 



as long as F u 3> a S> (see Fig. ||). This is true for dif- 
ferent types of distributions of the random increments 
{??(£, t)}, too. Except for the transients, the behavior of 
the model is also independent of the initial distribution 
of the F Cj o(«)' s on the interval [0, Our findings do 
not depend on the choice of boundary conditions as well. 
Periodic and open boundaries give similar power spectra. 
Even variations in the dynamic rules of the present model 
as realized, for example, in the Feder-Feder model |2j| do 
not alter our results. This points towards a generic be- 
havior strongly supporting the view in Refs. |ll|,[l6| that 
nonlinearity (here, activation/deactivation of sites with 
evolving thresholds) and dissipation are among the rele- 
vant features for generating 1// Q noise. 
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FIG. 3. log 10 -log 10 plot of the power spectrum of X{i) 
in our model with open boundaries for different cr's and 
(3 = 0.18, d = 2,N = 2500, F c = l,v = 0.1. The solid line 
with exponent 1.0 is drawn for reference. For d = 1 as well as 
for larger system sizes, we obtain similar results. 

The explanation of our results is the following: Due 
to the absence of critical behavior, only small avalanches 
occur in our model for dissipative /3's. Consequently, 
one can think of S(f) in a first approximation as the 
superposition of local power spectra S(f,i), 

N 

(8) 

j=0 

The local signal X(t,i) generating the respective S(f,i) 
is just the avalanche signal at site i, i.e., the number of 
topplings of this site during an avalanche at time t. This 
means that the sum over i of the X(t, i) is just X(t). 

We have investigated the local power spectra and we 
find indeed that Eq. (@) is a good approximation. This 
underlines especially that there is no dependence on the 
dimension of the system. It turns out that the S(f, i) 
are almost independent of i and that they show a 1//" 
behavior themselves (see Fig. |]). 

This result can be understood to a certain extent by 
mapping our model to a model introduced by Kaulakys 
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and Meskauskas p4| . In order to do so, we have to ne- 
glect all interactions between different sites. This means 
considering the limit (3 = or, equivalently, N = 1. Ad- 
ditionally, the random walk of the threshold is no longer 
confined by two reflecting boundaries. A parabolic po- 
tential centered around F c and characterized by the re- 
laxation rate 7 is used instead, 



F r 



F r + AF r 



AF C , T = (1 - 7 )Ai 71 c ,r-i + oti(t), 



(9) 



(10) 



with AF Ct o = 0. Since we consider N = 1, the time signal 
X(t) simplifies considerably: g.j = 1 for all j and, due to 
the uniform driving, the tj'a are given by 



Fr 



tj — tj-l 
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Fc) + *T)U), 
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(12) 



with F Ci i = and to = 0. Hence, X(t) is already deter- 



mined by the series of At, 



tj-i which evolve ac- 



cording to a random walk in a parabolic potential. This 
corresponds to one particle moving in a closed contour 
with the period of the drift of the particle around the 
contour fluctuating about the average value ^. The t, 's 
are then the transit times measured at a certain point. 
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FIG. 4. log 10 -log 10 plot of the power spectrum of X(t,i) 
in our model with open boundaries for different sites i and 
(3 = 0.22, a = 0.04, d = 2, TV = 2500, F c = l,w = 0.1. For 
bulk sites, a = 1.2 which is exactly the same as for X(t). For 
boundary sites, a — 1.1. 

This is indeed the model introduced in Ref. [gjj . 
Kaulakys and Meskauskas (KM) were able to compute 
analytically the power spectral density and obtained a 
power law with a = 1. This behavior can be found in any 
desirably wide range of frequencies for a sufficiently small 
7. The crucial point is that the law of large numbers is 
not valid for linin^oo — X)j l =o ^0 ( see ^ ne Appendix) . 
Their results explain at least qualitatively our findings in 
the limit (3 = 0. However, the exponent a differs. In our 



model, we find a = 0.91 ± 0.02 for = 0. This is due 
to an important difference between the KM model and 
our model. As already noted, 7 has to be small to gen- 
erate 1//" noise. Moreover, the asymptotic distribution 

2 

of the At's is a Gaussian with mean and variance 
This means that there is a non neglectable probability of 
negative At's. Consequently, tj + i can be smaller than 
tj implying a causality backwards in time. Hence, the 
KM model is somewhat ill-defined. In our model, nega- 
tive At's are not possible since the threshold has to be 
larger than or equal to zero. This also implies that the 
1// behavior cannot be extended to any desirable wide 
range of frequencies. 

In conclusion, we have shown that uniform driven SOC 
models are generally not able to generate l/f a noise 
without further modifications. The essential ingredient 
that has to be added is a local memory. This proves 
that spatially extended systems based on activation- 
deactivation processes do not have to be point driven 
to produce 1// Q noise. 

In the present model, the local memory is realized in 
the easiest possible way by a random walk of the thresh- 
old. We showed that the dynamics of the threshold is 
equivalent to the KM model of transit times under cer- 
tain assumptions in the limit (3 — * 0. This means that our 
model can be considered as a physically reasonable gen- 
eralization of the KM model to systems with a threshold 
- even if they are spatially extended. As a consequence, 
the present model combines the idea that 1// Q noise may 
result from the clustering of the signal pulses Q with the 
view that an activation-deactivation process and dissipa- 
tion are the main features relevant for the description of 
1// Q noise |ll^|l6| . The robustness of our results strongly 
supports these views. 

As an experimental realization of our model, we sug- 
gest a stick-slip system with a Markovian threshold evolu- 
tion. Finally, we would like to point out that the present 
model is similar to the coupled "integrate-and-fire" os- 
cillators studied in the context of neuronal networks and 
biology. Work is in progress to investigate these connec- 
tions further and might lead to an explanation for the 
occurence of 1// Q noise in cortical neurons || . 

J. Davidsen would like to thank the Land Schleswig- 
Holstcin, Germany, for financial support. 



APPENDIX: COMPUTATION OF THE POWER 
SPECTRAL DENSITY 

Consider a signal X(t) = ■ 6(t ~ tj) as in the KM 



model. Define Atj = tj 



spectral density [see Eq. (|l])] 



It follows for the power 
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With 7 = limT^oo ™(j 
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jmin + 1), this leads to 



5(/) = /(j]ex P l27r/9 ^^ i 
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where (...) denotes the average over the ensemble and 
over j. Hence, all we need is the probability distribution 
<f of the l/qJ2l~o ^ or tne model with an av- 
erage period At, it was shown that \& is a Gaussian with 
mean At and variance |- for / > j 3 / 2 /na pi] . Hence, 
l/ty y^ fc Aifc obeys the same distribution as Aij and does 



not depend on g. 
Eq. Q: 



This can be used to further simplify 



5(/) = J^(exp i2 ^E 
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For small enough /, the summation can be replaced by 
an integral. (In the KM model, this is valid for / < 
2^/7/ 7r<J an d / "C (27rAi) -1 .) Changing variables from q 
to q' — qf and evaluating the integrals gives 



S(f) = M(0)/f. 



(A6) 



The fact that the probability distribution of ^ X]fe=o ^ k 
does not depend on q merely means that the law of large 
numbers is not valid due to the correlation between dif- 
ferent Atfc's. 

This is of course different in the case of independent ran- 
dom variables Aft, i.e. for 7 = 1. In this case the dis- 
tribution of i _Cfc=o does depend on q. To be more 

precise, the distribution is still a Gaussian with the same 

2 

mean as before but the variance becomes . In the limit 
q — > 00, the variance vanishes. A short calculation gives 
as expected 
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P + fo 

This is exactly what happens in the KM model for 
/ > 7 3/2 /7TCT [||. 

To summarize, the crucial ingredient to obtain 1// 
noise is the generation of strong enough correlations be- 
tween different At, such that 1/q XZfc ^* obeys the same 
distribution as Atk does. The specific form of ^ is not im- 
portant as long as *£(0) 7^ 0. In the KM model, the strong 
correlations are implemented by a random walk dynam- 
ics. In general, other stochastic mechanisms are capable 
of generating such correlations as well, e.g., shot noise in 
combination with fast relaxation giving rise to random 
flows of events without memory and Cauchy statistics 
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